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Abstract 

The mechanism underlying the divergence of perturbation theory is 
exposed. This is done through a detailed study of the violation of the 
hypothesis of the Dominated Convergence Theorem of Lebesgue using 
familiar techniques of Quantum Field Theory. That theorem governs 
the validity (or lack of it) of the formal manipulations done to gener- 
ate the perturbative series in the functional integral formalism. The 
aspects of the perturbative series that need to be modified to obtain 
a convergent series are presented. Useful tools for a practical imple- 
mentation of these modifications are developed. Some resummation 
methods are analyzed in the light of the above mentioned mechanism. 
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1 Introduction 



A typical quantity to analyze the nature of the perturbative expansion in 
Quantum Field Theory is the partition function 



The normalization factor 1/Zq is the partition function of the free field 
{Z — > 1 when A — > 0). The analysis of the perturbative expansion of any 
Green's function goes along similar lines as in the case of Z. In the example 
above we consider a scalar field theory for simplicity. 

The traditional argument for understanding the divergent nature of the 
perturbative expansion can be traced back to Dyson Although the form 
was different, the content of his argument is captured by the following state- 
ment: 

"If the perturbative series were to converge to the exact result, the func- 
tion being expanded would be analytic in A at A = 0. But the function (Z 
for example) is not analytic in A at that value. Therefore, as a function 
of A, the perturbative series is either divergent or converges to the wrong 
answer." 

Estimations of the large order behavior of the coefficients of the per- 
turbative series showed that the first possibility is the one actually realized 
|^. |3|]. That Z, as a function of A, is not analytic at A = 0, can be guessed 
by simply noting that if in its functional integral representation (Eq. ([[])) 
we make the real part of A negative, the integral diverges. In fact, there is 
a branch cut in the first Riemann sheet that can be chosen to lie along the 
negative real axis, extending from A = — oo to A = Q, ||. 

The above argument is very powerful and extends to the perturbative 
series of almost all other nontrivial field theories. It has also motivated 
a series of very important calculations of the large order behavior of the 
perturbative coefficients , general analysis on the structure of field theories 




(1) 



with 
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pi], as well as improvements over perturbative computations of different 
physical quantities ||]. 

For all its power, it is fair to say that the argument, as almost any other 
reductio ad absurdum type of argument, fails to point towards a solution of 
the problem of divergence. It is only through the indirect formalism of Borel 
transforms that questions of recovery of the full theory from its perturbative 
series can be discussed || ]g, ||| . 

In this paper an alternative way of understanding the divergent nature 
of the perturbative series is presented. This way of understanding the prob- 
lem complements the traditional argument briefly described above, hopefully 
illuminating aspects that the traditional approach leaves obscure. In par- 
ticular, as we will see, the arguments in this paper point directly towards 
the aspects of the perturbative series that need to be modified to achieve a 
convergent series. It is hoped that the way of understanding the problem 
presented here will help to provide new insights into the urgent problem of 
extracting non-perturbative information out of Quantum Field Theories. 

In section || we develop our analysis of the divergence of perturbation 
theory. In Sec. |3| we point out the ingredients that, according to the anal- 
ysis of Sec. ||, a modification of perturbation theory would need to achieve 
convergence. We also present a remarkable formula ([34|) that allows us to 
implement such modifications in terms of Gaussian integrals, paving the way 
to the application of this convergent modified perturbative series to Quan- 
tum Field Theories. The proof of the properties of the function (|64[) is done 
in Appendix 1. In section || we analyze recent work on the convergence of 
various optimized expansions [12-19] in terms of the ideas presented here. 
In Sec. [| we summarize our results and mention directions of the work cur- 
rently in preparation. Finally, in Appendix 2, we apply the ideas of this 
paper in a simple but illuminating example for which we actually develop a 
convergent series by modifying the aspects of the perturbative series pointed 
out by our analysis as the source of divergence. 
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2 Lebesgue's Dominated Convergence Theorem and 
Perturbation Theory 



2.1 The wrong step in perturbation theory 

Although the notation will not always be explicit, we work in an Euclidean 
space of dimension smaller than 4 and in a finite volume. 

Let's remember how the perturbative series is generated in the functional 
integral formalism for a quantity like Z: 




The final sum is in practice truncated at some finite order N. The func- 
tional integrals that give the contribution of every order n are calculated 
using Wick's theorem and Feynman's diagram techniques with the corre- 
sponding renormalization. 

We see then that the generation of the perturbative series in the func- 
tional integral formalism is a two step process. First @ the integrand is 
expanded in powers of the coupling constant, and then the sum is inter- 
changed with the integral 

It will be convenient to have a simpler example in which the arguments 
of this paper become very transparent. Consider the simple integral 

2 (A) = _L f°° dxe-^+tO (6) 

A/7T J-oo 

1 In this paper we will often use the familiar word "integrand" to refer to e~ s or any 
functional inside the functional integration symbol. It would be more precise to preserve 
this word for e _Slnt in the measure defined by the free field. The terminology used here 
is, however, common practice in the Quantum Field Theory literature and also helps to 
emphasize the similarities with the intuitive finite dimensional case presented below. 
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and its corresponding perturbative expansion 

,„ i r. ? (-1)" i\ A" . x i ... 

= ' f f J^fWV (8) 

oo 

EE J2(-l) n CnX n . (9) 

n=0 

This simple integral has been used many times in the past as a paradigmatic 
example of the divergence of perturbation theory ||. It is then specially 
suited for a comparison between the traditional arguments and the ones 
presented in this paper. 

Again we see the two step process to generate the perturbative series. 
First the integrand is expanded in powers of A (0) and then the sum is 
interchanged with the integral @ . In this simple example the perturbative 
coefficients can be calculated exactly for arbitrary n. In the large n limit 
they become: 

y/2 

c n ~ — — (n — 1)! when n — > oo. (10) 

With such factorial behavior, the series diverges for all A different from 
zero as is well known. On the other hand the function z(X), as defined in 
Eq. (]6|), gives a well defined positive real number for every positive real A. 
Therefore one or both of the two steps done to generate the perturbative 
series must be wrong. 

Similarly, in the functional integral case normalized with respect to the 
free field (|l]), Z is a well defined number while its perturbative series diverges. 
One or both of the two steps must be wrong. 

The first step, the expansion of the integrand in powers of A, is clearly 
correct. As the integrand (not the integral!) is analytic in A for every finite 
A, the expansion merely corresponds to a Taylor series. The second step, 
the interchange of sum and integral, must therefore be the wrong one. 

The next obvious step is then to recall the theorems that govern the 
interchange between sums and integrals, to understand in detail why this 
is wrong in our case. The most powerful theorem in this respect is the 
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well known theorem of Dominated Convergence of Lebesgue. In a simplified 
version, enough for our purposes, it says the following: 



Let /at be a sequence of integrable functions that converge point- 
wisely to a function / 

►/ as TV^oo (11) 

and bounded in absolute value by a positive integrable function 
h (dominated) 

\f N \ <h, VN. (12) 

Then, it is true that 

lim I f N = I hm f N = [ f. (13) 



As a special case, if the convergence (11) is uniform and the measure of 
integration is finite, then the interchange is also valid. It should be empha- 
sized that Lebesgue's theorem follows from the axioms of abstract measure 
theory. Therefore if the problem under consideration involves a well defined 
measure, as is the case for the Quantum Field Theories considered here [19], 
the theorem holds. 

In our case we can write formally^], 

In [*(*)] = £ tf- (-J dV 4 ) V /dMi^tf+imV] (14 ) 

71 = 

for the functional integral case, and 

N 



n=0 

for the simple integral example. 

One important aspect of the dominated convergence theorem approach 
to analyze the divergence of perturbation theory is that it focuses on the 

2 See previous footnote. 
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integrands, objects relatively simple to analyze. On the contrary, the analyt- 
icity approach briefly described in the introduction focuses on the integrals, 
that are much more difficult to analyze. So, before we try to understand 
what aspects of the dominated convergence theorem fail in our case, let's see 
the "phenomena" (the integrand) for the intuitive simple example. In figure 




12 3 4 

X 



Figure 1: Exact integrand, zero , second and fourth perturbative approxi- 
mations. A = 1. 

|l|, the exact integrand, together with some perturbative approximations, are 
displayed. We can appreciate the way in which the successive approxima- 
tions behave. For small x, and up to some critical value that we call x Cj at, 
where the subindex c stands for critical while the subindex N indicates that 
this value changes with the order, the perturbative integrands approximate 
very well the exact integrand. Even more, x Cj at grows with N. But for x 
bigger that x C) n a "bump" begins to emerge. The height of these bumps, 
as we will see in detail shortly, grows factorially with the order, while the 
width remains approximately constant. So, the larger the order in pertur- 
bation theory, the larger the region in which the perturbative integrands 
approximate very well the exact integrand, but the stronger the upcoming 
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deviation. As we will see shortly, it is precisely this deviation that is respon- 
sible for the divergence of the perturbative series and the famous factorial 
growth. We will also see that an exactly analogous phenomena happens in 
the functional integral case and is again the responsible for the divergence 
of the perturbative series. 

Returning to the problem of understanding the aspect of the dominated 
convergence theorem that fail in the perturbative series we will now show 
that the sequence of integrands of Eq. (14) and Eq. (^) converge respectively 
to the exact integrands 



and 



F = — e~J dx i* (^r+fm^j-t J d u *r (i 6 ) 



f = 1 e -(*M* 4 ) (17) 
v vr 



but not in a dominated way. This is, there is no positive integrable function 
h satisfying the property of Eq. <K%). 



2.2 Failure of domination in the simple example 



That the sequence of integrands of Eq. (|14|) and Eq. ( |15| ) converge respec- 
tively to the exact integrands ( fl6| ) and flT7| ) is obvious, since, as mentioned 
before, for finite A they are analytic functions of A and so their Taylor ex- 
pansion converge (at least for finite field strength) . To see the failure of the 
domination hypothesis it is convenient to analyze the "shape" of every term 
of /at. Namely, for the field theory case, 

Cn [*(*)] = Y^T&l ^ I V/ d M^) 2 +X^j (1 8 ) 

while for the simple integrand 



c n (x) = ^-^(±Yx^e-* 2 . (19) 



w /tt n\ \4. 

In this section we analyze the failure of the domination hypothesis for the 
simple example @ because, as it turns out, it is remarkably similar to the 
Quantum Field Theory example analyzed in the next section. In Fig. 
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we can inspect the functions 03(2) and c^(x) for A = 1 corresponding to 
the simple integrand case that we analyze first. The maximum of c n (x) is 
reached at 

x max = ±(2n) 1 / 2 . (20) 
There, for large n, the function takes the value 

c n (* max ) = ^(-l)>-l)! A". (21) 

On the other hand, the width remains constant as n increases as can be seen 
by a Gaussian approximation around the maximum x max = (2n) 1 / 2 : 

ra ^72 (-!)> " X ) ! ex P [~ 2 ( x ~ ( 2n ) 1/2 ) 2 } A " ( 22 ) 
The integration of this Gaussian approximation gives, for large n 

dxc n (x)^^(-lT(n-l)\\ n (23) 

in accordance with Eq. ( [TCP if we take into account the factor of 2 coming 
from the two maxima ±(2n) 1 / 2 . 

The mechanism of convergence of the /jv's to / becomes clear now. The 
/at's are made out of a pure Gaussian (the "free" term) plus "bumps" (the 
perturbative corrections) oscillating in sign (see Fig. ^). The maxima of 
these bumps grows factorially with the order while their width remain ap- 
proximately constant (more specifically, the Gaussian approximation around 
the maxima (Eq. (|22|)), that becomes exact when the order goes to infinity, 
has a variance independent of the order). For fixed N, and for x smaller 
than a certain value, the bumps delicately "almost" cancel each other, leav- 
ing only a small remnant that modifies the free integrand into the interacting 
one. However, for x larger than that value, the last bump begins to emerge 
and, being the last, does not have the next one to get cancelled (in the 
N — > 00 limit, there is no last bump and the convergence is achieved for 
every x). Consequently, after a certain value x c ,n, the function /at deviates 
strongly from / and is governed by the uncanceled -/V th bump only, of height 
proportional to (N — 1)1 and finite variance. This is so because, since the 
height of the bump grows factorially with the order, for N large enough the 
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Figure 2: Exact integrand, and fourth perturbative approximation together 
with the third, and fourth terms. A = 1. 

last bump is far greater than all the previous ones, so it remains almost 
completely uncanceled. Since the variance of the bumps is independent of 
the order, this means that for every finite order, there is a region of finite 
measure in which the perturbative integrand is of the order of the height of 
the last bump. In figure || we can see how the function c±{x) is left almost 
completely uncanceled by c%{x) and dominates the deviation of f± from /. 

That x Cj tv (the value of \x\ up to which the perturbative integrand very 
accurately approximate the exact one) grows with N, going to infinity when 
N —* oo, is a simple consequence of Taylor's theorem applied to the analytic 
function e _Aa;4//4 . 

The above analysis makes clear the failure of domination of the sequence 
of Eq. (|j~5|) towards / (Eq. (|l7|)). Indeed, any positive function h(x) with 
the property 

IM*)| < h(x) , ViV (24) 

fails to be integrable, since it has to "cover" the bump, whose area grows 
factorially with N. So, although the sequence of /at(x)'s converges to f(x), 
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the convergence is not dominated as we wanted to show. 

Eq. (|23|), together with the above comments, indicate that the same 
reason for which the sequence of integrands ([H]) fails to be dominated, is 
the one that produces the factorial growth in the perturbative series. 

In the Field Theory case, although we can not rely on figures like || and 
U to guide our intuition, as we will show now, the analogy with the simple 
integral example is so close that the interpretation is equally transparent. 



2.3 Failure of domination in Quantum Field Theory 



For Quantum Field Theory, as for the simple example analyzed above, it 
is convenient to consider every term c n [</>(x)] (Eq. 18) of the perturbative 
approximation /jv (Eq. |l4|) of the exact integrand (Eq. [l^) , 



1 (-1) 



o n\ 



J d d 4i(9 M </>) 2 + imV]+nln[(A/4) / d d x <t> A 



(25) 



where we have written the n th power of the interaction in exponential form. 
The mathematical analysis below follows closely the discussions in chapter 
38 Ref. ||]. Although the problem treated there is different from the one 
treated here, many techniques used in || can be directly borrowed here. 

For n large enough, the analysis of its "shape" reduces to the familiar 
procedure of finding its maxima, as in the case of the simple integrand. The 
equation determining the maxima of c n [</>(x)] is the equation that minimizes 
the exponent, and can be thought of as the equation of motion of the effective 
action 

"A 
4 



s[<f>] 



d d x 



\ (<V) 2 + \m 2 4> 2 



nln 



d d x 



(26) 



which is 



V 2 6 + m 2 c 



An 



I d d xc, 



(27) 



Making the change of variables 



m 



An 



1/2 



ip(mx) 



m 



d/2-l 



4n 



fd d uip 4 (u) 



1/2 



ip{mx), (28) 
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<p satisfies the equation 

- V 2 ip(u) + tp(u) - ip 3 (u) = , u = mx. (29) 

This equation corresponds to the instanton equation of the negative mass 
Xcj) theory. The analysis of their solutions can be found in many places. We 
are interested in the solutions with minimal, finite action. For these solu- 
tions, in the infinite volume limit, scaling arguments provide very interesting 
information. We mentioned at the beginning of section 2.1 that we work in a 
finite volume. However, if the volume is large enough, the following infinite 
volume arguments remain valid up to errors that go to zero exponentially 
fast when the volume goes to infinity. 

Since the solution 4> max (x) (the subindex "max" indicates that, in func- 
tional space, c n \<j>{x)] reaches its maximum at (f>j agx (x), this should not be 
confused with the fact that the the action (|26| ) reaches its minimum there) 
is a minimum of the action (|26|), given an arbitrary constant a, S [a(/> max (x)] 
should have a minimum at a = 1 || [|. This implies the equation 

J d d x (d^ max ) 2 + m 2 J d d x<f> 2 max - 4n = 0. (30) 

Similarly, S [(pm^iax)] should also have a minima at a = 1, implying, 
(2-d) 



d 



J d d x {d^ max f -m 2 J d d x^ ax + 2n = 0. (31) 



Solving the system of equations (Bot) and (31) we obtain, 



d d x (^0 max ) 2 = nd (32) 



m 2 d d x< ax = n(4-d). (33) 



d d x 



2n, (34) 



From which we conclude in particular 

2 (d/^max) 2 + ^m 2 <t> 2 max 

independent of the dimension. The relations (|32|) and ( |33| ) can be explicitly 
checked in the case d = 1 (Quantum Mechanics), in which the solutions to 
eq.(^) are known analytically. They are 

^ {t) = O coshMi-io)] (35) 
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giving 



J dt(^) 2 = n (36) 
J dt^lf = 3n. (37) 



m 2 



Since ip(u), introduced in eq.(p8|) and satisfying eq. fl29|) , is dimension- 
less (remember that u = mx is also dimensionless) , and the corresponding 
0max(aO has finite action, the quantity 

A = 1 J d d n^(n) (38) 

is a finite, pure number greater than zero ||. For the Quantum Mechanical 
case mentioned above, A = 4/3. For the cases d > 1, A is not explicitly 
known but, as just said, it must be a finite, positive, pure number. With 
the definition (p8|), Eq. fl28|) becomes, 

0max(^) = m^ 2 ' 1 (^\ ' <p(mx) (39) 

Since (p(mx) satisfies the n-independent Eq. (|29|), we conclude that the field 
strength of ^ max grows with the square root of the order n. 

Equation ([34]) , together with the definition (|38|) and the relation 
allow us to write an expression for the action ( [26] ) at <p = </>max> 



5* f^maxl = 2n - n In 



A m d " 4 2 



.4 

The value of c n [<j>{x)\ at <p = (/> max then becomes, for large n, 



(40) 



n 



Cn[^ m ^)]^Y ^ in ~ iy \^A~) ■ (41) 

With the change of variables 

<t>{x) = 0max(x) + m^ 2 -^ q (mi), (42) 
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the Gaussian approximation of c n [(f)] around max is, 

. i^-^m"- 



e 2 

_ i 
e 2 



§/ d d Wld d W 2 ^q(«l)[(-V 2 (1 +l-3 ¥ . 2 (u 1 ))(5(«l-«2)]^ q («2) . (43) 
|J d d M 1 d d M 2 <A q («l)[(l/ J 4) V 3 («i)^ 3 (M 2 )]0q(«2) (44) 



where u = mx and <p(u), solution of Eq.(29), is related to </> max through 
Eq.(p9|). This Gaussian approximation becomes exact in the limit n — > 00. 
The second derivative operator, that we call D, is then, 

D D\ oc& \ + D non _ local 

(45) 

with 

Aocal = -V 2 + 1 - 3if 2 (46) 

and ^ 

Awn-local = -j^v ><v \, with < u\v >= f 3 (u) (47) 



and A given in Eq.(|3£ 

The operator D\ oca \ is well known (see for example |J). It has d eigen- 
vectors \0fj_ > with zero eigenvalues given by 

<«|0m>=^(«)- («) 

These vectors are also zero-eigenvectors of -D, as can be seen by noting that 
\v > is orthogonal to them, 

< V \Q >= 0. (49) 



They reflect the translation invariance of the action (26). 

-f local is also known to have one and only one negative eigenvector. The 
proof of this fact given in Appendix 38 of reference Q, that uses Sobolev 
inequalities, can be repeated line by line to prove that, on the contrary, D 
is a positive semi-definite operator, 

D > (50) 
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in the operator sense. 

Projecting out the d-dimensional eigenspace of eigenvalue zero, the re- 
sulting operator, that we call D', is positive definite. 

D' = A'ocal + Amu-local > (51) 

This equation explicitly states that the projection over the strictly positive 
eigenvectors modifies only Diocal- The non-local part, as we saw, is a projec- 
tor orthogonal to the zero modes and is therefore not modified under that 
operation. 

Equations ( |50| ) and (|5l|) suggest that the operator D, with the corre- 
sponding renormalization for d > 1, generate a well defined Gaussian mea- 
sure in a finite volume (remember d < 4). In fact, the determinant of D [ ocal 
was calculated many times in the past ||, and a generalization of a Quan- 
tum Mechanical argument of ref. [10] indicates that this is all we need to 
compute the determinant of D'. The argument goes as follows, 



A'ocal + -j\v >< V\ 



Det [D 1 ] = Det 

= Det [A'ocal] (l + J< v \ D 'iocJ v >) ■ (52) 
Since ip(u) is orthogonal to d^tp(u) (the zero modes of D and Mocai), 

A'ocal = Aocal^ = "V- (53) 

The last equality follows from the definition of -Diocal in Eq.(|46"|) and the 



equation (29) satisfied by tp. Inverting A'ocal' anc ^ remembering the defini- 
tion of \v > and A in Eqs. ( f47[) and (|38|), we obtain 



< v\D^\v >= -2A. (54) 



Replacing this result in Eq.(52|), we arrive at the result 

Det [A] = -Det [D[ ocal ] . (55) 

As already mentioned, Aocal nas one an< ^ orn y one negative eigenvector, con- 
sequently its determinant is negative. Eq. (J55|) indicates then that Det [D'\ 
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is positive, as it should be according to ([5l]). The effect of the nonlocal part 
is to change the sign of the determinant of the local part. 

The preceding equations allow us to integrate the Gaussian approxima- 
tion of c n \<p(u)\ given in Eqs. ( fi3|]4i| ). Using the method of collective co- 
ordinates to project out the zero modes, the Jacobian of the corresponding 
change of variables is, at leading order in 1/n, 



J 



d 

n 

n=i 



\2 jd 



fj, rmax J 



1/2 



(56) 



where no sum over n is implied. 

It can be shown that the solutions of Eq. (|2~7| ) corresponding to minimal 
action are spherically symmetric S, then (|5"6|) can be written as 



J 



! 2 d d x 



-id/2 



(57) 



where now, sum over \i is implied. Using Eq. (32) we then find 

J = n d / 2 . (58) 
With this expression, the functional integral of c n [<p(u)] can be written as 

{-IT, .jAm^" 



1 

Zq 



IT / [<*0] C n [</>] 



2vr 



(n- 1)! 



(Vol m d ) 



n 



d/2 



A 



-Det 



D' 



local 



Da 



-1/2 



(59) 
(60) 



where-Dg = —V 2 + 1. The factors in the line ( |59j ) correspond to the value 
of c n [eft] at (ft maiX up to the normalization 1/Zq as can be seen in Eq. (|4l|). 
The factor "Vol" arises after the integration over the flat coordinates cor- 
responding to the center of max - The n d l 2 comes from the Jacobian of 
the change of variables as mentioned before. The factor m d arises after the 
rescaling of the fields that makes them dimensionless in both c n [(f)] and Zq. 
This happens because there are d more integration variables in Zq due to 
the integration over the collective coordinated in the numerator. Finally, the 
factor (—Det [-D( ocal ]) ' is the result of the integration over the coordinates 
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1 /2 

orthogonal to the zero modes of D, while (Det [-Do]) is the dimensionless 
normalization factor (the mass dimension of both, the numerator and the 
denominator, was already taking care of in the term m d ). The minus sign 
is due to the non-local part of D that, as proved above, simply changes the 
sign of the determinant of the local part, making it positive. 

In the case d = 1, -Bet [D[ oc J D ] = 1/12 §, [U|, and Eqs. f>j||60|) 



(with A = 4/3 as already mentioned) become identical to the corresponding 
result of Ref. [ 10 1 if we take into account the different normalization here 



and a factor of 2 that is taken care of by remembering that the sign of the 
solution of Eq. ( |27| ) is undetermined, therefore both, positive and negative 
solutions contribute equally to the functional integral. 

For d = 2 or 3, the formal expression (^J6^) needs of course to be 
renormalized. All the arguments in this section remain valid for the theory 
with a Pauli-Villars regularization ||. The action @ becomes 



S [0] = J d d x 




+ ^ 4 + i<W(A) 



(61) 



The modification of the kinetic part of the action affects both the equa- 
tion (|27]) and the scaling arguments, but by an amount that becomes small 
like A -2 when the ultra-violet cut-off A becomes large. 

As shown in Ref. , although the counterterm increases with the cut-off, 
since it is proportional to at least one power of A, taking the small A limit 
before the large cut-off limit justifies to ignore it in the equation ( p7|) and the 
scaling arguments. On the other hand it contributes to the result ( ^S|j60| ) an 
amount that exactly cancels the divergence in the Det [-Di oca i] making the 
final expression finite as it should be. 



In the large n limit, where the Gaussian approximation (|43|, 44) becomes 
exact, the expression (^J6^) gives the large order behavior of the pertur- 
bative series of Z (up to the factor of 2 mentioned above) without any as- 
sumption about the analytic structure in A [pLCtl ■ A completely analogous 
procedure would give the large order behavior of any Green's function. 

Eqs. (J3^) , (PI) , (I431J441) and (l59| , |60|) , allow us to draw an accurate pic- 
ture of the mechanism underlying the lack of domination (in the sense of 
the Lebesgue's theorem) of the convergence of the sequence of perturbative 
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integrands ( |l4| ) towards (|l6|) , and consequently of the mechanism underlying 
the divergence of the perturbative series. In fact, perhaps not surprisingly 
given the similarity of their large order behavior, this picture is very similar 
to the one described in the previous section for the simple integral example. 

In a finite volume, there is a region of finite measure in field space in 
which the perturbative approximation f^[(j)(x)] of Eq. ( |T4| ) approximate the 
exact integrand (|l^) with an error smaller than a given prescribed number. 
This region grows with N, becoming the full field space in the N — > oo limit. 
As in the simple example, this is a consequence of Taylor's theorem applied 
to the (analytic) integrand (16). 

The problem is that, for any finite N, outside that region the approxi- 
mate integrand f]\i[(p(x)] strongly deviates from the exact one. This can be 
seen by noting that the maxima of every term of /jv grow factorially with 
the order. Therefore, for large enough N, the last term is far greater than 
the previous ones at its maxima. Even more, as shown above, the Gaussian 
approximation around that maxima (that becomes exact for N — > oo) de- 
fines a measure that does not go to zero as N — > oo (in fact, it is independent 
of N ( 43,^4])). This means that for every finite N, there is a region of finite 
measure in field space (and this measure does not go to zero as N — > oo) in 
which the deviation between the perturbative integrand and the exact one 
is of the order of the maxima of the last term of _fjv, i.e., of the order of 
(N — 1)!. No integrable functional can therefore satisfy the property ( |l2| ) of 
the Lebesgue's theorem. 

That is the mechanism that makes the sequence of perturbative inte- 
grands, although convergent to the exact one, non-dominated in the sense 
of Lebesgue's theorem. That is therefore the mechanism that makes the 
sequence of integrals (i.e., the perturbative series) divergent. In fact, as 
Eqs. (59,60) show, the famous factorial behavior of the large order coeffi- 
cients of the perturbative series is a consequence, after integration, of the 
above mechanism. 
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3 Steps Towards a Convergent Series 



It was mentioned in the introduction that the analysis of the divergence of 
perturbation theory presented in this paper would point directly towards 
the aspects of the perturbative series that need to be modified in order to 
generate a convergent series. This is the topic of the present section. 

In the previous section we analyzed perturbation theory from the point 
of view of the Dominated Convergence Theorem. We have detected the 
precise way in which the convergence of the sequence of perturbative in- 
tegrands to the exact one takes place, and the way this convergence fails 
to be dominated. We have learned that for any finite order N, the field 
space naturally divides into two regions. In the first one, that grows with 
the order, eventually becoming the full field space (in the N — > oo limit), 
the perturbative integrands very accurately approximate the exact one. In 
the other one, however, the deviation between the perturbative and exact 
integrands is so strong, that the sequence of integrals diverge. 

It is then clear that if we could somehow modify the integrands, order by 
order, in the region where they deviate from the exact one, while preserving 
them as they are in the other region, then, with a "proper" modification, 
such modified sequence of integrands would converge in a dominated way. 
According to the Dominated Convergence Theorem, their integrals would 
then converge to the exact integral, achieving the desired goal of a conver- 
gent modified perturbation theory. 

Let us call VLjq to the region of field space in which the N th perturbative 
integrand approximate with a given prescribed error the exact integrand 
(|l6|). The characteristic function, Ch(£lN, {(fr(x)}), of that region is equal 
to 1 for field configurations belonging to it, and zero otherwise: 



One possible realization of the above strategy of modifying the integrands 




for {4>(x)} G Q, N 
for {(t){x)} g" Q.N 



(62) 



( |l4]) in the "bad" region of field space is to make them zero there. We would 
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have 

i N (—l} n /A r \ n 

f' N M*)] = Yo E ~nT e ' S ° (l J Ch( ^' (63) 

According to the analysis of the previous section, choosing £In appro- 
priately, the sequence of f' N [4>{x)\ will converge dominatedly, and the cor- 
responding interchange between sum and integral will now be allowed. A 
rigorous proof of this is left for a paper currently in preparation. For the 
purposes of the present argument, it is sufficient to rely on the analysis of 
the previous section to assume its validity. Also, in the next section we will 
analyze, along the general ideas of this paper, some resummation schemes 
for which rigorous proofs of convergence have recently been given [12-19]. As 
that analysis will show, these methods strongly rely on the general notions 
underlying the formula (|65|). Their convergence supports, then, the validity 
of the dominated nature of the convergence of (^) towards ( 1C ) . 

An urgent issue, however, is the practical applicability of the above strat- 
egy. To implement it, we need a functional representation of the character- 
istic function ( |62] ) (or an approximation to it) that only involves Gaussian 
and polynomial functionals. In the same way in which a functional represen- 
tation of the Dirac delta function allow us to perform functional integrals 
with constraints, the Fadeev- Popov quantization of Gauge theories being the 
most famous example, a functional representation of the characteristic func- 
tion ( |62| ) would allow us to functionally integrate only the desired region of 
functional space. Since, basically, the functionals we know how to integrate 
reduce to Gaussians multiplied by polynomials, the desired representation of 
the characteristic function should only involve those functionals. Conversely, 
if it only involves those functionals, all the sophisticated machinery devel- 
oped for perturbation theory (including all the perturbative renormalization 
methods) would automatically be applicable. With this in mind, consider 
the following function, 

W(M,u) = e- Mu J2 { -^- (64) 



7! 



where M is a positive integer. Note that W(M, u) arises from 1 = e e 



Mu„+Mu 
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by expanding the second exponential up to order M. W(M, u) has the fol- 
lowing remarkable properties, 



1. W(M,u) — ► 1 when M — > oo for < u < 1. The convergence is 
uniform, with the error going to zero as 

R(M,U) < e M(ln«-(u-l))_l = ( g 5) 

V ^ ~ v 7 ^? 1 - u + 1/M v ; 

2. W(M, u) — ► when M — > oo for 1 < it. The convergence is also 
uniform , with an error of the form 

W(M,u) < e ^Cinu-(tt-i)) < (66) 

As we see, the exponent corresponds to the same function in both cases. 
For u > 0, this function is always negative except at its maxima, at u = 1, 
where it is 0. Therefore the convergence is in both cases exponentially fast 
in M, with the exponent becoming more and more negative, for a fixed M, 
when u differs more and more from 1. The proof of properties 1 and 2 is in 
the Appendix 1. 

If we replace u by a positive definite quadratic form < 4>\D\4> > /Cjv, 
then the insertion of Eq. (|34|) into the functional integral would effectively 
cut off the region of integration < 4>\D\<fi > > CV 

= ^E^Jta / [d „e-(5,„ t )^(M,lil^ 
^o^, n! M~*coJ \ C N 

(68) 

Cat is a constant that changes with the order of the expansion in A, 
increasing with iV but in such a way that in the region < (f>\D\<p >< CV 
the difference between the perturbative and the exact integrands is smaller 
than a given prescribed error. Since the convergence of W is uniform ac- 
cording to properties 1 and 2, with errors given in Eqs. ( |65|) and (|66|), the 
corresponding interchange between the sum in Eq. (p8|) and the functional 
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integral is justified. The fact that u becomes a quadratic form implies that 
the resulting integrands are Gaussians multiplied by monomials, therefore 
the familiar Feynman diagram techniques can be used to integrate them. It 
also implies that no new loops appear and the sum in j from ( |64|) becomes 
an algebraic problem. A typical functional integral to compute has the form 

J [d 0] e -/d^[|(a^) 2 +|m^^( W /c7 JV )] d d X(j A n d d x m y 

(69) 

as can be seen by replacing the definition ( |64|) into ( |68| ) with u =< <p\D\<fi > 
/Cat- 

Note that at any given order in A, it is not necessary in principle to go 
to infinity in M. That would amount to replace the perturbative integrands 
by zero in the region < (f>\D\<p > > Cn, realizing the strategy mentioned 
before. But since the convergence in W is uniform, a finite, large enough M 
(depending on the order in the expansion in the coupling constant), would 
suffice to tame the behavior of the perturbative integrands and transform 
them into a dominated convergent sequence. In fact, as we will see, many 



methods of improvement of perturbation theory use effectively formula ( 64 ) 
without sending M — > oo for any given finite order in perturbation theory. 
In any case, as already mentioned, that limit is in principle computable, 
since it does not involves new loops. Work in this direction is in progress. 



The convergence of the sequence (68) towards Z(\) may be thought, at 



first sight, to be in conflict with our well established knowledge about the 
non-analyticity of this function at A = 0. In fact, Eq. (68) seems to be a 



power series in A (the powers of A coming from the powers of Si n t), there- 
fore, if convergent, that power series would define a function of A analytic 
at A = 0. It must be recognized, however, that the validity of Lebesgue's 
Dominated Convergence Theorem is completely independent of any analyt- 
icity consideration. Therefore, if its hypothesis are satisfied, its conclusions 
must be valid. This being said, the question of how does the convergence 



of (68) fits with the non-anlayticity of Z{\) deserves an answer. To begin 
with, even at finite order in A, the function ( |68|) is not necessarily analytic 
at A = despite its analytic appearance. This is because the constant Cn 
may have an implicit nonanalytic dependence on A. In Appendix 2 this is 
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actually the case in the context of a simple example to which the present 
ideas are applied. But the mechanism that ultimately introduces the proper 
non-analyticity in A is the limit process N —* oo. Given a non- analytic 
function like Z{\) one can always construct a sequence of analytic functions 
that converge to it. Satisfying the hypothesis of the Dominated Conver- 
gence Theorem is a way of achieving that, avoiding all the complicated and 
model dependent issues of non-analyticity. Note that the validity of these 
hypothesis for a given sequence of integrands can be checked independently 
of any analyticity consideration. 

In the Appendix 2 we prove the convergence of the general strategy dis- 
cussed here for the simple integral example analyzed in section (|2.2|). For 
that case, making u = (x/x C) n) 2 , the function W(M, u) becomes in the 
limit the characteristic function of the interval \x\ < x Cj jv. We use this 
to explicitly compute the non-analytic function z(X) (Eq. @) calculating 
only Gaussian integrals. We also show explicitly how a non-analytic de- 
pendence of x Cj tv on A naturally arises just by demanding the validity of 
the Lebesgue's hypothesis and how the N — > oo limit process captures the 
full non-analyticity of z(X). The same method also works for the "negative 
mass case", where the Borel resummation method fails. In figure ||, we can 
appreciate the convergence of W towards the characteristic function of the 
interval \x\ < x Cj tv for x Cj tv = 1 for two different values of M. 



4 Improvement methods of perturbation theory. 

The analysis of the mechanism of divergence of the perturbative series pre- 
sented in this paper, together with the formula ( |64| ) and its properties, offer 
a large range of possibilities to construct a convergent series. In the previous 
section we have shown how that formula can be used to effectively cut-off 
the region of field space where the strong deviation between perturbative 
and exact integrands take place. But as we will see, this is only one possi- 
ble way, among many, to use the formula (^34|) to transform the sequence of 
perturbative integrands into a dominated one. 

Another example of its possible use is the so called "optimized delta 



expansion" pO, ET[. In a series of papers [ 11 1- [ 14 ] , it was proved that such an 
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Figure 3: Function W(M, x, x C) jv) with x c ^ = 1 for M = 3 (segmented line) 
and M = 60 (continuous line). The convergence towards the characteristic 
function of the interval \x\ < x c n is apparent. 



expansion converges for the partition function of the anharmonic oscillator in 
finite Euclidean time. The problem of convergence in the infinite Euclidean 
time (or zero temperature) limit for the free energy or any connected Green's 
function is still under investigation, as well as its extension to Quantum Field 
Theories [|l^, 15 1. The method was proved to generate a convergent series 
for the energy eigenvalues [pi 17 1, although such studies make heavy use of 



analyticity properties valid specifically in the models studied. In these works, 
it was realized that many methods of improvement of perturbations theory, 



such as the order dependent mappings of references |22], 23 1 , posses the same 
general structure as the linear delta expansion. A considerable amount of 
work has been dedicated to investigate the virtues and limitations of the 
method and extensions of it fl 



It is not the place here to give a detailed analysis of these methods. But 
we would like to briefly indicate how they can be understood in terms of 
the ideas presented here. In what follows, our analysis is restricted to d = 1 
(Quantum Mechanics) where rigorous results about the convergence of the 
methods considered here are available. 

Let's consider the case of the anharmonic oscillator. Its action is given in 
Eq. (Eh for d = 1. The idea of the method is to replace it by an interpolating 



24 



action 



Ss 



dt 



m 2 + - — a 
Zm 



a 

— ( 

m 



(70) 



Clearly, the dependence on the parameter a in 5,5 is lost when 5 = 1. For 
that value, the action ( |70"| ) reduces to (g). However, if we expand up to a 
finite order in 5 and then make 5 = 1, the result still depends on a. The 
idea is to tune a, order by order in the expansion in 5, so that the result 
is a convergent series. It was shown in the references mentioned above that 
the methods works if a is tuned properly. For example, in reference fi"2[| , 
the asymptotic scaling a ~ jV 2 / 3 was used to prove the convergence of the 
method for the partition function at finite Euclidean time. 

], a was tuned according 



It is interesting to note that originally [11 



to heuristic prescriptions such as the "principle of minimal sensitivity" [gl]] 
(at any given order in 5, choose a so that the result is insensitive to small 
changes in it), or the criterion of "fastest apparent convergence" (the value 
of a at which the next order in delta vanishes). But later fl3| , [H| , it was 
realized that the best strategy was simply to leave a undetermined, find an 
expression for the error (that obviously depends on a), and then choose a 
so that the error goes to zero when the order in 5 goes to infinity. It is clear 
that a structural understanding of the convergence of the method can help 
to construct the necessary generalizations to overcome the difficulties associ- 
ated with the convergence in the infinite volume limit for connected Green's 
functions, as well as the extensions to general Quantum Field Theories. 

To understand the "optimized delta expansion" in terms of the ideas 
presented in this paper, let us expand the functional integral corresponding 
to the action (70) in powers of 5 up to a finite order N, and make 5 = 1 as 
the method indicates, 



Z(m,X,a,N) 



1 

N 

£■ 

,n.=0 



-Jd t [i(d t 0) 2 + I( m 2 + ^)^ 



4 



Xa 
4m 



(71) 



The general analysis of the mechanism of divergence of perturbation 
theory of section indicates that if the function ((TTj) generates a convergent 
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series with a scaling properly with N, then, barring miraculous coincidences, 
the corresponding integrands should converge dominatedly (or, even better, 
uniformly) towards the exact integrand ( |I4| ) . We want to obtain a qualitative 
understanding on how this method achieves that. 

Expanding the binomial and making some elementary changes of vari- 
ables in the indices of summation, we obtain the expression 



Z{m,X,a,N) = — [ [dcf>] e ~ I dt ^ {dt(l ' )2+ ^ m 



2 , \a ' 




(72) 



This equation already shows some of the distinctive characteristics of the 
method. As we see, the i th power of the interacting action in the expansion 
of e~ Slnt up to order N, is multiplied by 



W (N-i) = e-™^f* £±(£[<lA . (73) 




Note that W (N) corresponds to the function W(M, u) with M = N (N 
is the order in the expansion of e~ Slnt ), and the variable u replaced by the 
quadratic form ((A/4m) / (f> 2 ) /Cat, where CV = N/a. Making for example 
a ~ N 2 / 3 as in Ref. (where it was proved that with such scaling the 
method generates a convergent series), we see then, that, according to the 
previous section, W (N) is an approximation of the theta function in the 
region of field space characterized by 

^ fdx <P 2 < N 1 / 3 . (74) 
4m J 



Equation (72), however, shows that the mechanism used to achieve dom- 
inated convergence can not be reduced to a simple insertion of the function 
W(M, u) with M = N and u = ((A/4m) / 4> 2 ) /C N . That would be the case 
if all the powers of the expansion of e~ Int up to order were multiplied 
by W (N). But equation fl7"2"D shows that the i th power of the interacting 



action is in fact multiplied by W (N — i). 
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At this point it is convenient to pause for a moment in our study of the 
"optimized delta expansion" to give some useful definitions. 

Let us call passive mechanisms (to achieve dominated, or uniform con- 
vergence of a sequence of integrands to the exact one) to those that can be 
reduced to the product of the -/V th perturbative integrand and the charac- 
teristic function of a region £1^ of field space for some sequence {^at}- 

Passive methods use only information that is already available in the 
perturbative integrands, they just get rid of the "noise" inherent to pertur- 
bation theory. Because of that, in addition to define a convergent series, 
they can also be very useful to study perturbation theory itself. The func- 
tion W(N, u), with u replaced by a properly selected quadratic operator, was 
specially designed to make passive methods practical. In a sense, section || 
is a discussion of passive methods. 

Active mechanisms are those that are not passive, as defined above. 

What kind of mechanism is the one underlying the "optimized delta 
expansion" method? 

A trivial generalization of the proof, in the previous section, of the con- 
vergence of W(M, u) towards the theta function for u > 0, shows that the 
function 

W(M, u, i) = e~ Mu £ (75) 

n=0 l - 

also converges towards the theta function for u > in the limit 

M -» oo, i fixed. (76) 

In this sense, the "optimized delta expansion" method does have passive 
aspects. As Eqs. (fT2j f73|) show, it amounts to multiplying the i th power of 
the expansion up to order N of e Slnt by W with u = ((A/4m) / (j) 2 ) /Cn and 
Cn = N/a. Since this function converges to the characteristic function of 
the region characterized by Eq. (|74|), this means that the first i terms, of 
the expansion up to order N of e Slnt are effectively multiplied by the same 
function (an approximate characteristic function) for i <S N . Therefore, 
the first i terms, with i <C N, use only the information available in the 
perturbative series to converge to the exact integrand. 
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What about the other terms?, the ones characterized by i < N? Sur- 
prisingly, these terms produce a convergence of the corresponding integrands 
towards the exact one that is faster than possible with only passive compo- 
nents! 

It is not the place here to study this aspect in detail, so, let us simply 
see this "faster than passive" convergence for the simple integral example. 

Applied to the "mass-less" version of the integral (||) , the optimized delta 
expansion method was proved to generate a rapidly convergent sequence in 
Ref. iPHl . That is, the sequence given by 

lN = j2 t^i r dx e -«w - 2 ( ^4 _ a{N) x 2 )n (77) 

n=0 n - J -°° 4 

was proved to converge to 

/oo 
dx e~ Xx /4 (78) 
-oo 

when a(N) ~ y/N with an error that goes to zero at the very fast rate of 
Rn < C'A^ 1 / 4 e _0 ' 663Ar when N — ► oo. C is a numerical constant. 

We are interested in understanding whether the corresponding conver- 
gence of the integrands is faster than passive. For our qualitative purposes, 
it is enough to observe, in figure ||], the convergence towards the exact inte- 
grand 

/cxa(x) = e"^ 4 (79) 
of both, the perturbative integrand 

w - ± ^ m 

n=0 

and the optimized delta expansion integrand 

N (_-i\n \ 

/ode = £ M^e-^W x \-x± - a{N) x 2 ) n (81) 

with a(N) ~ y/N, for N = 4. 

We can see how accurate the convergence of I d e (x) is, already at this 
low order. In particular, when the perturbative integrand begins to diverge, 
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Figure 4: In the main plot, the superiority of the convergence of the fourth 
order optimized delta expansion (ode) with respect to the same order per- 
turbative approximation is evident. In the sub-graph, the difference between 
the ode and the exact integrand is plotted. Note the difference in the scales 
of the y axis of the main and sub graph. 



29 



Iode(x) continues to approximate remarkably well the exact integrand. In 
the sub-figure, we can appreciate the difference between I eX a(x) and / ot j e (x). 
Note the difference in the y axis scale of graph and sub-graph. 

It is then clear that the optimized delta expansion method, with its 
subtle combination of passive and active components, manages to generate 
a sequence of integrands that (uniformly) converges towards the exact one 
at a rate that far exceeds the possibilities within a purely passive method. 

We stop here the qualitative discussion of the optimized delta expansion 
method expressing two general lessons: 

1. Any method of improvement of the perturbative series in a given quan- 
tum theory, where a functional integral representation of the quantity 
under study exists, must rely, at the level of the integrands, on an im- 
provement over the pointwise convergence of the Taylor series in the 
coupling constants of e~ s . 

2. The problem of finding a convergent series reduces to the problem of 
finding a dominatedly convergent sequence of integrands towards e~ s . 

This second simple statement, not only provides a guide to the construction 
of convergent schemes, but also emphasizes the fact that, in principle, a 
dominatedly convergent sequence of integrands, do not have to have any 
relation whatsoever with the corresponding Taylor expansion. In order to 
be able to use the usual techniques of Quantum Field Theory, it is reasonable 
to restrict the search for a convergent scheme to a sequence of integrands of 
the general form: 

N 

f N = e~ So a n SL fn(< <P\D n \4> >) (82) 

n=0 

where the functional f n of the quadratic form < (p\D n \<j) >, should take 
care of the non-dominated convergence that is bound to appear with only 
powers of the interacting action. The function W of section [3|, with its 
possible generalizations, is an ideal candidate for this purpose. But the 
selection of the coefficients a n amounts to a pure problem in optimization of 
the convergence of the integrands — no a priori connection with any Taylor 
series is necessary. 
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5 Conclusions 



In this paper we have exposed the mechanism, at the level of the integrands, 
that make the perturbative expansion of a functional integral divergent. We 
have seen in detail how the sequence of integrands violate the domination 
hypothesis of the Lebesgue's Dominated Convergent Theorem. That theo- 
rem, as is well known, establishes the conditions under which it is allowed 
to interchange an integration and a limit, in particular the one that takes 
place in the generation of perturbation series. 

It was shown that at any finite order in perturbation theory, the field 
space divides into two regions. One, that grows with the order, in which the 
perturbative integrands very accurately approximate the exact integrand. 
In the other however, a strong deviation takes place. It was shown that the 
behavior in this second region violates the hypothesis of Lebesgue's theorem, 
and, consequently, generates the divergence of perturbation theory. The 
famous factorial growth of the large order coefficients of the perturbative 
series was shown to be an effect, after integration, of the very mechanism 
that violates the hypothesis of the theorem. 

All of the above was done explicitly without relying in the particular an- 
alytic properties of the models studied. It is therefore natural to assume that 
similar mechanisms of violation of the Lebesgue's hypothesis are present in 
any other Quantum Field Theory, although for just renormalizable theories 
other mechanisms are responsible for renormalons. Studies in this direction 
are in progress. 

The mechanism of divergence presented here points towards a simple 
way to achieve a convergent series: integrate only in the "good" region of 
field space. Since this region grows with the order, becoming in the limit 
the whole field space, integrating in a correspondingly increasing region we 
would obtain a convergent series. A step forward towards a practical imple- 
mentation of this program was made with the construction of the function 
W (|64|). This function allows us to introduce a Gaussian representation of 
the characteristic function of regions of field space, very much like the im- 
position of constraints in the functional integral was allowed by a functional 
representation of the Dirac's delta function. A rigorous proof of the conver- 
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gence of this practical implementation of the above mentioned strategy is in 
progress. In Appendix 2 it was applied to a simple integral example. 

Finally, a qualitative analysis of the optimized delta expansion method of 
improvement of perturbation theory in terms of the ideas of this paper was 
done. Some general properties of improvement methods, useful to generates 
new schemes, as well as to understand and improve old ones, have been 
established. 
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Appendix 1 

In this appendix we will prove the two properties of formula (|64|). 

Because for u > all the terms of the sum defining W(M, u) are positive, 
we have trivially W(M, u) > 0. On the other hand, since in the Taylor 
expansion of e Mu all the terms are positive, we have J2^?=o (Mu) n jn\ < e Mu . 
Therefore W(M,u) < 1. So for every M and positive or zero u we have, 

0<W(M,u)<l. (83) 

Consider first the case < u < 1. 

l-W(M,u) =e- Mu Y ^^=R(M,u), (84) 

n=M+l 

we will prove that R(M, u) — > when M — > oo. 
Changing variables to j = n — M, we get 

R{M, u) = e-^^t { M U y^ (85) 
< e -M„(M^gJ^t (g6) 



-Mu(Mu) 



M 



£ e M\ l-u+l/M - (87) 



But M M /Ml -> e M /V2^M for large M, so 
R(M,u) < e ^(ta«-(«-i))_A 



v^M l-u + l/M' 
The exponent is negative in the region < u < 1 since, being both, In u and 
(n — 1) negative there, | lnu| > \u — 1| in this region. Therefore 

R(M, u)->0, when M oo (89) 
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in the region < u < 1 and property 1 is proved with an exponentially fast 
convergence. 

In the region u > 1, we have 

W(M,u) = e-Mu^iMuT^^M^Ml {90) 

n=0 n=0 
< e M(ln«-(w-l))_ /g-Q 

The first inequality is valid because u > 1 and the second because 
e M > Xm=o M n /re!. The exponent is again negative. For u > 1, both, lnu 
and (u — 1) are positive, but now | lnu| < \u— 1|. So property 2 is also valid 
with an exponentially fast convergence. 

For u = 1 all we know is that W is bounded by Eq.([33"|). That is all we 
need. Numerics suggest W(M, 1) — > 1/2 when M — > oo. 

This finishes our proof. 



Appendix 2 

In this appendix we apply the strategy discussed in section || to generate a 
series convergent to the function z(X) (Eq. (||)). This is done using the func- 
tion W of Eq. (p4]) and computing exclusively Gaussian integrals, therefore 
we restrict ourselves to using only those techniques that are also available 
in Quantum Field Theory. 

As mentioned in section ||, the simplest possible modification of the per- 
turbative integrand ( |l5| ) that would transform the corresponding sequence 
into a dominated one, amounts to keep them as they are for \x\ < x c ^ and 
replacing them by zero for \x\ > x Cj tv- That is, 

jl _ J ^ Z^n=0 ^ 4 X J e IOT l X l < (92) 

i for \x\ > x c ^n 

In fact, choosing x C) n so as to properly avoid the region in which the devia- 
tion takes place, the sequence of f' N converges uniformly towards the exact 
integrand (|l7|) as we will show shortly. Consequently, the corresponding 
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sequence of integrals 



dxf' N = tt-^/ ■ dx^M^ > 4 ) e " X ( 93 ) 

-*c,JV n=0 " : V ' 



iV 



n 



will converge to the desired integral 

z (A) = tt- 1 / 2 f° dxe-^ 2 ^ 4 ). (95) 



In (93) the change in the limits of integration from ±oo to ±x C) n is just due 
to the definition of f' N in Eq. ( |9~2| ) . The interchange between sum and integral 
in (|94[) is now allowed because in the region [—x c ^n,x c ,n] we have uniform 
convergence (this is a stronger condition than dominated convergence) . The 
resulting integrals are not Gaussian due to the finite limits of integration. 
We will show how they can be calculated using only Gaussian integrals. 

A trivial way to achieve convergence of the sequence of integrals of the 
f' N of Eq. (^|) towards (|95|) amounts to keep x Cj n equal to a finite constant 
''a" independent of N, while taking the limit N — > oo. In this limit, Eq. (94) 



becomes identical to 7r -1 / 2 J** dxe~( x ' + 4 X ), since for finite a the Taylor 
series of the integrands converge uniformly. Therefore, as already said, the 
interchange between sum and integral is legal. Finally, taking the limit 
a —* oo, we would obtain the desired convergence towards z(X). 

However, better use can be made of the information available in f' N for 
finite N. For example, for every finite N, we can choose x c ^ so that 

\f N ( x )- f( x )\<^L for \x\<x C)N , (96) 

with €t,n going to zero as N — ► oo. Then, since we have 

\f' N (x) - f(x)\ < e -H"+K") = ^ for \x\>x c>N , (97) 



the fjtf(x) will uniformly converge towards the exact integrand f(x) if ( p6| ) 
is consistent with x c ^ — ► oo when N — ► oo. Indeed, if this happens, we 
would have 

POO 

1/ (f{x)- f N (x))dx\ <e TjN + e c , N ^0 when N -> oo. (98) 
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The term €t,n comes trivially from (^), while e C) N comes from Eq.(|97|) and 
the inequality 



e 



< x2+ -^)dx < e -«" + *<") = e Cj7V , (99) 



valid for x Cj at > 1. 

Applying Taylor's theorem to the function e _Ax ' 4//4 one can easily show 
that the condition (]96T) is satisfied if 



X C ,N 



(iV + 1)! 



e T ,N ^4\^ +1 )l 1/(4(7V+5/4)) 

aJ 



(100) 



Note that the non-analytic dependence of x c ^ on A arises automatically 
from the imposition of Eq. ( |96| ) to satisfy the hypothesis of the Lebesgue's 
theorem. 

Remember that the only condition on tT N to achieve convergence of 
the sequence of integrals is to go to zero when N — > oo consistently with 
x c,N — ► oo in that limit. Choosing for example 

e T ,N = e-* Nl/ \ (101) 

we obtain asymptotically, 

x CjiV ^(4iV/eA) 1/4 . (102) 
This implies (through Eq. (|97|)), 

ec N ^ e -meX?"-N/e t (1Q3) 



Equations ( |101| ) and (103) show the exponential rate at which the conver- 
gence of the sequence of integrals take place. 

Clearly the form (|101| ) for ex,N is not unique, not even the most efficient 
one, but enough to achieve convergence. 

In the table (|[) one can appreciate the numerical convergence for A = 
4/10. 

Up to now we have proved that the general strategy of section [3| does, in 
fact, generate a convergent sequence towards .z(A). However, the resulting 
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Table 1: Integration over the small field configurations only produces a 
convergent series. In the last column the improvement over the perturbative 
values can be appreciated. 



Order 


Exact value (A = 4/10) 


Conv. series 


Pert, series 


2 


0.837043 


0.803160 


0.848839 


4 


0.837043 


0.830264 


0.854087 


6 


0.837043 


0.835516 


0.901897 


8 


0.837043 


0.836667 


1.316407 


20 


0.837043 


0.837044 


2.33755 10 8 



integrals in (^4|) are not Gaussians, making the applicability of the method 
in Quantum Field Theory dubious, to say the least. We will show now that 
the integrals of Eq (p4|) can be computed, using Eq. ( |6"4| ) with u = (x/x Cj jv) 2 , 
calculating only Gaussian integrals. The steps involved are 

/Xc,N o roc _ 

x r e~ x dx = / x r e~ x lim W(M, x, x c , N )dx (104) 



co 2 



lim / x r e~ x W(M,x,x cN )dx (105) 



3C 

M 



^to n] \ x In) J-™ 

(106) 



The two properties of W validate both equality (104) and (because of the 



uniformity of the convergence in W) dlQ5| ). In the last line (^0|) we just 



make explicit the meaning of (|105| ). So it is clear that these two properties 
are enough to prove the validity of ( |106j ) , where only Gaussian integrals are 
present. But it is a good exercise to find a direct proof of it in the case at 
hand, where everything can be computed exactly. We do this next. 

For r odd the integrals vanish, so let's consider the case r even, that is, 
r = 2t, for any integer t. 

On the one hand we have 
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where the necessary interchange between sum and integral to arrive to the 

2 

result is allowed due to the uniform convergence of the Taylor series of e~ x 
in the finite segment [—x Ct N, x c ,n]- 
On the other hand, 



lim > 



1 



M 



M^oo 



n=0 

M 



n\ \ x 



c,N , 



-(l+M/xl N y x 2(n+t) da 



lim J2—T(n + t+ 1/2) 



n=0 
M 



lim > 



n 



l x 



L c,N 
M 



t+1/2 



1 + 



x: 



N 



M 



c,N 



t+1/2 oo 



n=0 



fc=0 



A:! 



-T(n + t + k + 1/2) 



M 



(108) 
(109) 
(110) 



(x c ,N) 2t+1 



-1) 



(x c ,n) 



2k 



k=0 



k\ (k + t + 1/2) 



(fc + t + 1/2) T (n + t + fc + 1/2) 



n=0 



n! 



(Ill) 



In line fll09j) we have used the equation 



r>*\ h: r(n + l/2) 



j; e 



,n+l/2 



(112) 



in line (110) we have expanded the last term of (109) in powers of x 2 cN /M 
and carried out some cancellations, and finally in (111) we have interchanged 
the M — > oo limit with the infinite sum in k. 

Comparing (107) and ( |111| ), we see that the validity of Eq. ( p.06| ) depends 
on the validity of the equation 

(k + t + 1/2) JH, T(n + k + t + 1/2) , . N 

lim V ,uZ,L Y , = 1 v integers k, t > 113 

n=0 

That this identity holds for every integer t and k can be seen by considering 
the following analytic function of the complex variable z: 



0{z) 



(1/z 



M 



lim 

M-coo M\ 



■ ^) \ - *_\ 
i/z) 2^ T 



n=0 



T(n+l/z) 
'(n + l)! ' 



(114) 



If the identities (113) hold, this function must be identically 1, since for 
1/zj = j + 1/2 with j integer it reduces to them, and for ever increasing j, 
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we obtain a sequence accumulating at z = on which the function should 
be 1. 

Conversely we will prove that 0(z) is indeed identically 1 as an analytic 



function of z, proving in consequence the identities (113) for arbitrary t and 
k. Consider the sequence 1/zj = j + 1 for j integer. This sequence also 
accumulates at z = 0, and for all its points we have 



0(1/(7 + 1)) 



l im O' + l) Y> T(n + j + l) 
M (i+i) r(n + l)! 



M^oo M<i +1 ) 

r (J + 1 ) 
hm — , . ■ 4 

(j + 1) 



= lim 



n=0 

M 

Ln=0 

(JTT) +0 < Mr 



M^oo 



(115) 
(116) 
(117) 
1 (118) 



Therefore O(z) = 1 for all z. This finishes the direct proof of Eq. (106). 

As was mentioned before, Eq. ( |1Q0|) , that was derived independently of 
any analyticity consideration, and only with the purpose of satisfying the 
hypothesis of Lebesgue's theorem, introduces a non-analyticity in the se- 
quence of integrals of f' N even for finite N. But even for the case where x c ^ 
is fixed to a constant a, discussed before, in which the limit N — > oo is taken 
first, and then a is sent to infinity, and therefore the sequence is made out of 
truly analytic functions, the convergence towards z(X) is perfectly compat- 
ible with analyticity considerations. The functions 7r -1 / 2 J** dxe~( x + * x ) 
(the result of the N — > oo limit), are clearly analytic in A. But they con- 
verge to (in fact they define!) the nonanalytic function z(X) when a — » oo. 
The limit of an infinite sequence of analytic functions does not have to be 
analytic. 

Another important issue is that the same method works also for the 
"negative mass case", where the Borel resummation method fails. Indeed, 
from the discussion of this section it must be obvious that, with a proper 
scaling of x Ci n, the f' N 's with negative quadratic part of the exponent also 
converges uniformly towards e^ x2 ~^ x ^ for x in [—X(^Ni x c,n]- Therefore, the 
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sequence of integrals is also convergent. 
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